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By treating CP-violating interaction of the electroweak bubble wall as a pertur- 
bative term, chiral charge flux through the bubble wall is estimated. It is found 
that the absolute value of the flux Fq has a sharp peak at uiq ~ a ~ T with 
Fq/(uT s ) ~ 10~ 3 (Ql — Qr) AO. Here mo is the fermion mass, 1/a is the wall thickness, 
T is the temperature at which the bubbles are growing, u is the wall velocity, Ql^r) is 
the chiral charge of the relevant left(right)-handed fermion and AO is the measure of 



CP violation. 
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1. Introduction 



Provided that the electroweak phase transition is first order with the bubble 
formation 1 ^ 2 ), a fermion gets a complex mass through the transition in the presence 
of CP violation: m(z) = mji(z) + imj(z) = mo(/(z) +ig(z)), where z is the coordinate 
perpendicular to the wall, f(z) gives the profile of the wall, mo is the fermion mass at 
the temperature in question, and the imaginary term g(z) yields CP violation necessary 
for baryogenesis 3 ) . 

In a previous paper, we have examined fermion scattering off the CP- violating bub- 
ble wall by treating g(z) as a perturbative term, and have derived the CP- violating effect 
AR = Rft^i — Rji—,.L AO in the kink-type background of f(z) = (1 + tanh(az))/2 
(distorted wave Born approximation — DWBA) 4 -*^. Here Rr^^Rr^i) is the reflec- 
tion coefficient for the right-handed fermion (anti-fermion) incident from the symmetric 
phase region and reflected as the left-handed one, and AO is the measure of CP- violation. 
| AR | is found to decrease exponentially as a /mo decreases, where 1/a is the thickness 
of the bubble wall.^ This means that, given 1/a, the heavier fermion whose Compton 
wave length is less than the wall thickness gives the smaller contribution to | AR |. In 
other words, a large Yukawa coupling does not necessarily cause a large effect of CP 
violation in electroweak baryogenesis. 

In this paper, we analyze the chiral charge flux Fq oc AO just in front of the bubble 
wall in the symmetric phase region. Fq is given by integrating AR weighted by the 
fermion flux densities in Sec. 2. In Sec. 3 numerical estimates of Fq are given. Fq has 
a sharp peak at mo ~ a ~ T, where T is the temperature at which the bubbles are 
growing. At the peak we find Fq/(uT 3 ) ~ 1CT 3 (Ql — Qr)A6, where u is the wall 
velocity and Qur) ls ^ ne chiral charge of the relevant fermion. Several remarks on 
electroweak baryogenesis are given in Sec. 4. 



f Reference 4) is referred to as I hereafter. 

ff Because of some errors in our numerical calculations, figures and table in I are incorrect. 
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2. Chiral Charge Flux through Bubble Wall 

Suppose that a left(right)-handed fermion of the species i with chiral charge Q l L r R y 
which is conserved in the symmetric phase, is incident to the wall from the symmetric 
phase region. The change of the expectation value of the chiral charge in the symmetric 
phase region, which has been brought by the reflection and transmission of the fermions, 
is given by 

AQ? = Wr - Q l L )RU R + (-QI + QrWr^l 

+ (-Qi)(n^ L + iu R ) - {-Q R ){f s R ^ L + f^ R )] r Li 

+ [(Qi - Qr)R r ^ l + (-Q R + Qi)RU R 

+ (-Qr)(T^l + Rr^r) ~ (-Qim^L + T[^r)] fm- 



(2.1) 



The notations would be obvious: T's are the transmission coefficients and f[m^ is the 
left(right)-handed fermion flux density in the symmetric phase region, and /£ = f R has 
been taken into account. Making use of unitarity relations examined in I, we obtain 

AQ? = (Qi - Q R ) [(Rr-^l ~ RU R )f! + T^Xfh - f s Ri )}, (2.2) 

where ff = (1/2) (/^ + f Ri ) and is the transmission coefficient in the absence of 
CP violation. Note that 

AR = R S R ^ L - R S R ^ L (2.3) 

and Aff = ff^ — f Ri are quantities of 0(A6). In a similar way to (2.1), the change of 
the expectation value of the chiral charge brought by the transmission of the fermions 
incident from the broken phase region is 



aq^ = q1(tL l /L + Tju/jy + QWi^fu + Tidily 

+ (-Ql)(T b R ^ L f b Li + T^Lfm) + (-QR)(T b R ^Rf b Li + f h L -> R fm 
In addition to unitarity relations in the broken phase region, 

T L-,L + T L-,R + R b L^L + R L^R = 
Tr-^l + Tr^r + Rr^l + Rr->R = 1) 



(2.4) 



(2.5) 
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the following reciprocity relations hold 5 ); 



rpb I rpb rpS , rpS I 1 pS \ 

1 L^L + 1 R-^L - 1 R^L + 1 R~^R\- 1 ~ K R^L)i 



(2.6) 

Thanks to these relations, we obtain 



AQ* = -(Qi - Q R ) [(R S R ^ L - RU R )ff - * T^Uh - fm)} , (2.7) 

PL 

where pi is the longitudinal momentum of the fermions to the bubble wall. Putting all 
these together, the chiral charge flux just in front of the wall in the symmetric phase 
region is given by 

00 oo 

F l Q = Q \~ 2 ^ R J d PL J PTd PT [f!(PL,PT) ~ f!(p L ,PT)}&R(a/m ,p L ). (2.8) 



m 

Here 

ft = (p L /£?)(exp[ 7 (E - up L )/T] + I)" 1 , 

ft = (p L /E)(eMl(E - u^p\ - ml) /T] + l)" 1 , (2 ' 9) 
are the fermion flux densities in the symmetric and broken phase regions respectively, u 



is the wall velocity, 7 = 1/ \/l — u 2 , E = y p 2 L + p^, Pt being the transverse momentum 
of the fermions along the wall. AR is defined in (2.3). In (2.8) we have ignored the 
effect of CP violation on the fermion flux densities (the second terms in (2.2) and (2.7). 
See the discussions in Sec. 4). 
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3. Numerical Estimates of Chiral Charge Flux 



Let the fermion mass be 



m(az) 



m(x) = mo(f(x) + ig(x)) 



with 



x = az. 



(3.1) 



We take the kink-type profile of the wall with the thickness \ja: 



f(x) = (1 + tanhx)/2, 



(3.2) 



and apply DWBA to the imaginary part g(x) oc AO to calculate the CP- violating effect 
AR oc AO. The next order corrections would be O((A0) 3 ) since AR is a CP-odd 
quantity. 

The profile of the complex mass m(x) should be determined by solving the equations 
of motion of the gauge-Higgs system with the effective potential at the temperature 
considered. In the case of moderate first order phase transition, we expect the modulus 
of the mass to be the kink-type profile. Its phase may be determined by the CP-violating 
nature of the effective potential as well as the boundary condition in the broken phase 
region, which is restricted by experiments at zero temperature. To have the explicit form 
of g{x), we need to specify the model, to calculate the effective potential and to solve 
the equations of motion with an appropriate boundary condition which is consistent 
with experiments. This will be very important but difficult process to estimate the 
baryon asymmetry quantitatively. We do not execute this program here but assume 
two forms of g(x) to evaluate the chiral flux. Note that it is g'(x), not g(x) itself, that 
contributes to AR as a factor of the integrand of it.^ 

In the first case, g(x) is related to the real part f(x) in a simple manner: 

(a) g(x) = A0f n (x) with n = 2.^ 

In the second case, g'(x) has the same thickness with that of f(x) and tends to zero 
both in the symmetric and broken phase regions: 

f See Eqs.(2.25), (2.27) and (2.35) in I. 
ff AR = for n = 1 as pointed out in I. 
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(b) g'{x) = A#sech(x), the CP-violating range being equal to the wall thickness. 

These are free from any restriction from experiments, and give large | AR/A9 \. 

Fig.l shows AR/A8 for the case (a). For all forms of g(x) and g'(x) we have 
analyzed, | AR | as a function of a/ mo decreases as a/ mo decreases almost exponentially 
(thick walls). 

Fig. 2 illustrates u dependence of Fq/u for the case (a), which is found to be almost 
constant for u < 0.5. 

Fig.3 shows the contour plots of Fq/(uT 3 ) for u < 0.5 for the cases (a) and (b). The 
both plots share a common pattern. Namely, they have the crest lines along mo ~ a, 
and fall particularly rapidly for mo > a since AR itself decreases exponentially as stated 
above. The crest lines are sharply peaked at mo ~ a ~ T. There we have 

F«y( M T 3 )^ KT 3 (Qi-Q^AO. (3.3) 

— Figs.1,2,3 — 

4. Remarks on Electroweak Baryogenesis 

First we briefly review how the chiral charge flux transforms to the baryon number 
through the sphaleron process along the charge transport scenario 3 '*. We take the weak 
hypercharge Y as the conserved charge in the symmetric phase. Fermion distributions 
near the bubble wall are changed by the hypercharge flux, which are to be obtained 
by solving the Boltzmann equations. Suppose that the initial states are in chemical 
equilibrium with B = L = through the gauge and Yukawa interactions. New chemical 
equilibrium emerges by injection of the hypercharge flux. These situations are described 
by the chemical potentials fi. From the chemical equilibrium and initial conditions, the 
chemical potentials for the baryon number (hb) an d the weak hypercharge (/zy) are 
related as 

HB = -/iy/6. (4.1) 



G 



The hypercharge density ny is expressed by the chemical potential jiy as 



T2 i 1 + 2m ) (a o^ 

= y 4 V y, (4-2) 

where m is the number of Higgs doublets. The transformed baryon density p B in the 
symmetric phase region is obtained by integrating the detailed balance equation, 

Tb r B 6 

P-B = — ^r/^ = - — j^n B {z, t), (4.3) 

where T B is the sphaleron transition rate. Considering that n(z, t) depends only on the 
distance from the wall in the rest frame of it, we put n(z, t) = n(z — ut) = n(z'). Then, 

z/u oo 

PB = ~^7f2 J n B (z - ut)dt = ^—^-JL- J n Y {z')dz'. (4.4) 

-oo 

The last integral is estimated as 

oo 

- [ n Y (z')dz' ~ ^V, (4.5) 
u J u 





where r l is the transport time within which the scattered fermions are captured by the 
wall 3 ). 

If t % is given by D l /u with the diffusion constant D l it is estimated from 
L>*(quark)~ 1/T and ^(lepton)- (10 2 ~ 10 3 )/T. From Y B = na^T 4 ~ 10~ 6 T 4 with 
k — 0.1 — 1.0, we obtain 

? - 1 »- 7 >< £ 4- <«» 

where the entropy density is given by s = 27r 2 g*T 3 /45 with g* ~ 100. 

As we saw in the previous section, Fy strongly depends on a and thq. The effective 
potentials of the minimal standard model making use of the high-temperature expansion 
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suggest rather thick bubble walls (l/a = (20 ~ 40)/T) 8 ) 9 ) 10 ). We have for T ~ lOOGeV 
and u — 0.1, 

p B /s ~ 10~ 18 x AO (for top quark), 

(4.7) 

Pb/s ~ 10~ 13 x A6 (for r lepton). 

Some models, such as a two-Higgs-doublet model proposed by three of the authors 
(K.F., A.K. and K.T.) without resorting to the high-temperature expansion n ), predict 
thin walls with l/a ~ 1/T. Then we expect rather large F Y which leads to 

Pb/s = (10~ 9 ~ 10~ 10 ) x A6 (for top quark and r lepton). (4.8) 



Now some comments are in order. 

(1) D l could be estimated from the Boltzmann equations and is expected to be 
D l = P/3 where l % is the mean free path of the fermion. But there is a suggestion of an 
enhancement of r by a factor of O(10 3 ) due to the forward multiple scattering effects 3 ). 
If so, this leads to 

Pb/s = (10~ 6 ~ 10~ 7 ) x AO, (4.9) 

and would be able to explain the observed value. 

(2) As u — > 0, pb in (4.6) diverges because r*(= D % /u) — > oo. However, phys- 
ical considerations suggest that the limit u — > leads to a contradiction to the 
nonequilibrium condition. Some discussions were made that there is a cutoff veloc- 
ity u m in ~ 10~ 4 — 10~ 2 7 ) 9 ) and that pb — * as u — > at u < u m i n . 

(3) In the minimal standard model, the condition to suppress the sphaleron tran- 
sition after the electroweak phase transition leads to the Higgs mass bound m# < 50 
GeV in contradiction to the experimental results. Also in the model, the CP-violating 
effect emerges only in the KM matrix and is too small to explain the observed value of 
Pb/s. Extensions of the standard model by incorporating two Higgs doublets (including 
SUSY model) have the possibility to avoid the above upper bound of m#. 

In order to get F y quantitatively, however, we must know the spatial dependence 
of the CP-violating phase 9(x). We are trying to solve nonlinear equations of motion 
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for 6(x) of the two-Higgs-doublet model 5 ) 11 ) under some assumptions for the real part 
of the bubble wall. 

(4) A time- dependent 9(x,t) plays a role similar to the chemical potential and is 
called charge potential. It influences the flux densities /^m^ and causes baryogenesis 
through the sphaleron process (spontaneous baryogenesis scinario) 3 ). The function 
Aff = /£j — /L is given from 8(x,t) through the Boltzmann or diffusion equations. 
For thick walls, the second terms in (2.2) and (2.7), which have been ignored in this 
paper, would become effective. Some results under the classical and hydrodynamical 
treatment were reporte 

In conclusion, we need to obtain the profile of the bubble wall and, in particular, 
the functional form of 8(x,t) in order to estimate F Y quantitatively, which is the basic 
quantity in any scinario of electroweak baryogenesis. 
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Figure Cations 

Fig.l. AR/A8 in the case (a) vs. Pi/mo for several a/mo. 

Fig.2. Fy(uT 3 (Qi - Q* R )A6) vs. u at a/T = 10, 1, 0.1 in the case (a) with a/rriQ = 5. 

Fig.3. Contour plot of F^/(uT 3 (Q^ - (%)A0) at u < 0.5 ClS db function of the wall 
thickness 1/a and the fermion mass mo. (a) g(x) = A9 f 2 (x). (b) g'(x) = 
A9 sechx. 
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